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Abstract 

The Falicov-Kimball model consists of spinless electrons and classical particles 
(ions) on a lattice. The electrons hop between nearest neighbor sites while the ions 
do not. We consider the model with equal numbers of ions and electrons and with a 
large on-site attractive force between ions and electrons. For densities 1/4 and 1/5 
the ion configuration in the ground state had been proved to be periodic. We prove 
that for density 2/9 it is periodic as well. However, for densities between 1/4 and 
1/5 other than 2/9 we prove that the ion configuration in the ground state is not 
periodic. Instead there is phase separation. For densities in (1/5,2/9) the ground 
state ion configuration is a mixture of the density 1/5 and 2/9 ground state ion 
configurations. For the interval (2/9, 1/4) it is a mixture of the density 2/9 and 1/4 
ground states. 
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1 Introduction 



The spinless Falicov-Kimball model has two types of particles: spinless electrons and 
classical particles, which we refer to as ions. The particles are on a lattice with the 
restriction that there is at most one ion at each lattice site. The spinless electrons are 
fermions, so there is at most one electron per site. The electrons can hop between nearest 
neighbor sites, but the ions cannot. There is an on-site interaction between electrons and 
ions. The Hamiltonian is 

H= E ctc y -AUj2clc x V x (1) 

<x,y> x 

where c£ and c x are creation and annihilation operators for the electrons. V x is the 
occupation number for the ions, i.e., V x — 1 if there is a ion at x and V x = if there is 
not. The sum over < x, y > is over nearest neighbor bonds in the lattice. (The factor of 4 
in front of the U is included for latter convenience.) This paper is only concerned with the 
square lattice, although the model may be defined on any lattice. We will only consider 
the neutral model in which the number of electrons is equal to the number of ions, and 
the interaction parameter U will be large and positive. By a hole-particle transformation 
results for positive U imply results for negative U, but we will not bother to state them. 

A review of rigorous work on the Falicov-Kimball model may be found in 0. Here 
we mention only some of the work on the neutral model for large positive U. In one 
dimension it is expected that for large U the ground state of the neutral model with 
rational density is the periodic arrangement of the ions which is "most homogeneous." 
(There is an explicit algorithm for determining the most homogeneous configuration.) 
This was proved by Lemberger for U > U c where U c depends on the denominator of the 
rational density [0. In any number of dimensions the ground state for density 1/2 is 
the checkerboard configuration for all U > In two dimensions with large U the 

ground states for densities 1/3, 1/4 and 1/5 are known rigorously and are periodic || 
For densities between 1/4 and 1/2 there are partial results on the ground state 0, but 
there is no proof it is periodic for rational densities. Based on what is known in one and 
two dimensions and the methods used to obtain these results, it is natural to conjecture 
that in two dimensions the ground state for large U is periodic for rational densities. 

In this paper we prove that this conjecture is wrong for densities between 1/5 and 1/4 
other than 2/9. In this density range there is phase separation in the ground state. The 
phases involved are the ground states for densities 1/5, 2/9 and 1/4, which are shown in 
figure 1 . The ground state for densities between 1/5 and 2/9 is made up of large regions 
of density 1/5 and density 2/9 ground states with the relative areas chosen to yield the 
desired density. A similar statement holds for densities between 2/9 and 1/4. The precise 
theorem is as follows. 
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A B D 

Figure 1 : The ground states for densities 1/5, 2/9 and 1/4 (left to right). The 
heavy lines between the ions are only a guide for the eye. 



Theorem 1 There are positive constants Uq and c such that for U >Uq the following is 
true for L by L squares A with n ions and n electrons. Let p = n/L 2 be the common 
density. If the density p is 2/9, L is a multiple of 6 and we use periodic boundary 
conditions, then the ground state configurations of the ions are configuration B in figure 
1 and its translates. If the density p is between 1/5 and 2/9, then for every ground state 
configuration of the ions we can find a subset A of A which contains at most cU 8 L sites 
and is such that on each connected component of A\A the configuration agrees with either 
configuration A or B in figure 1 (up to a lattice symmetry). If the density p is between 
2/9 and 1/4, then the same statement is true with "A or B" replaced by "B or D." (The 
statements for densities other than 2/9 are true for any choice of boundary conditions.) 

Although U 8 is large, the key point is that the bound on the number of sites in Ao 
contains L while the number of sites in A is L 2 . So for large L, A is a tiny fraction of the 
total area. We can think of A as consisting of domain walls between large regions which 
contain one of the configurations shown in figure 1. The phase separation we find in this 
range of densities for the neutral model should not be confused with the phase separation 
discussed by Freericks and Falicov ||. Their argument applied to the non-neutral model. 

The expectation that the neutral model with rational density should have periodic 
ground states is based on the following intuition. The attraction between electrons and 
ions is large, so each electron spends most of its time at a site with an ion. Now consider 
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the kinetic energy of an electron. If nearby sites have ions, then the electrons at those 
sites will restrict the movement of the electron we are considering. So its kinetic energy 
is minimized by spreading out the ions as much as possible to maximize the space that 
each electron has to move in. However, Watson || emphasized that there will typically 
be a mismatch between the lattice and the natural ion configuration in the absence of 
a lattice. Thus the lattice structure can frustrate the exclusion principle's attempts to 
put the ions in the "most homogeneous" configuration. We should emphasize that this 
paper only covers a small interval of densities. An important open question is whether the 
phase separation we find here holds for most densities, or whether there are intervals in 
which the rational densities have periodic ground states which are the most homogeneous 
in some sense. 

We conclude the introduction by sketching the proof for densities in (1/5, 2/9). When 
U is large and the model is neutral, the ground state energy of a given ion configuration 
may be expanded in powers of 1/U. This yields an effective Hamiltonian for the ions. 
One can begin to study it by only keeping terms up to a certain order in 1/U. (Of course, 
to prove anything one must eventually consider all orders.) Watson |3j showed that when 
the density is between 1/5 and 1/4 the ground states of the fourth order Hamiltonian 
correspond to tilings of the plane by squares and diamonds in which the squares and 
diamonds have the dimensions of those found in figure 1 . Watson's result will play a 
crucial role in our proof. The vertices in such a tiling can be one of four types which we 
label A, B, C or D following Watson's notation. The four types are shown in figure 3 . 
Note that in figure 1 configurations A, B or D contain only vertices of type A, B or D, 
respectively. 

All configurations that correspond to a square-diamond tiling have the same energy 
through fourth order. To determine the ground states for densities in (1/5, 2/9) we must 
go to higher orders in the perturbation series. Following Watson's treatment of a similar 
model, we write the Hamiltonian as a function of the number of each type of vertex. At 
sixth order these square- diamond tilings still all have the same energy. At eighth order 
vertices of type A,B and C have the same energy but vertices of type D have higher 
energy. For densities between 1/5 and 2/9 there are square-diamond tilings which contain 
no vertices of type D. To determine the ground state among all these configurations we 
must go to tenth order. Here we find that the energy of a type C vertex is higher than 
that of a type B vertex. Thus for densities in (1/5,2/9) the ground state must be a 
square- diamond tiling with only type A and B vertices. (Note that the tiling with only 
type A vertices has density 1/5, while the tiling with only type B vertices has density 2/9. 
So one can obtain any density in (1/5, 2/9) by a suitable mix of type A and B vertices.) 
However, a type A vertex cannot be adjacent to a type B vertex. Thus we must separate 
the A and B vertices to minimize the energy. 
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2 Proof of phase separation 



To derive the perturbation theory it is convenient to change to "spin" variables for the 
ions. Let 2V X — S x + 1, so S x = 1 when there is an ion at x and S x = —1 when there is 
not. We then take 

<x,y> x 

This differs from the original Hamiltonian by a term proportional to J2x c t c x, but we 
will only consider problems in which the number of electrons is fixed, so such a term is 
constant. 

There are no interactions between the electrons, so the Hamiltonian is just the second 
quantized form of the single electron Hamiltonian T — 2U S. The operator T has matrix 
elements T xy with T xy = 1 if \x — y\ = 1 and T xy = otherwise. The operator S is 
diagonal with entries S x . The ground state energy for N electrons is the sum of the 
N lowest eigenvalues of T — 2U S. To find the ground state for a particular density of 
electrons and ions we must minimize this energy over all S with the desired ion density. 

Let H(S) be the ground state energy for the ion configuration S with the number of 
electrons equal to the number of ions. To expand H(S) in powers of 1/U, we begin by 
rewriting H(S) as in 0. If U > 2 then the number of negative eigenvalues of T — 2US is 
equal to the number of sites with S x = 1, i.e., the number of ions. Thus when the number 
of electrons equals the number of ions, we have 

H(S) = ]T A, = \[Tr(T - 2US) - Tr(\T - 2US\)] (1) 

where Aj are the eigenvalues of T—2U S. Now Tr(T) = 0, and if we keep the number of ions 
fixed then Tr(S) is a constant. So we might as well redefine H(S) = — ~Tr(|T — 2US\). 
Then we write this as 

H(S) = ~Tr(\T - 2US\) = -\tt{[{T - 2US) 2 }^) = -UTr{l + A) 3 (2) 

with 

A = -{2U)- l {TS + ST) + {2UY 2 T 2 

We have used the fact that S 2 = 1. 

Now we derive the perturbation theory by following the treatment by Gruber, Jedrze- 
jewski and Lemberger ||. A somewhat different derivation may be found in [[/J. If U is 
sufficiently large, then ||A|| < 1 and we may expand (1 + A) a in a power series in A. 
Since T xy is nonzero only if x and y are nearest neighbors, when we take the trace of each 
term we generate nearest neighbor walks that end where they start. Grouping together 
terms with the same power of U , we may write the result as 

oo 

H(S)= £ U- 2m+1 J2 h 2m , x S x (3) 

m=l X 
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X is summed over finite subsets of the lattice, and Sx = Tlxex &x- The coefficient h 2m ,x 
is nonzero only if there is a nearest neighbor walk with 2m steps which ends where it 
begins and visits each site in X. (It may visit sites outside of X as well.) The coefficients 
h2m,x are invariant under the usual lattice symmetries. There is a constant c such that 
for every site x 

E IWI<c m (4) 

X:xeX 

We need to compute this effective Hamiltonian through tenth order, i.e., through the 
m = 5 terms. There are a lot of terms at this order. Since we will consider relatively 
low densities, it will prove useful to go back to the occupation variables V x . (Recall that 
S x = 2V X — 1.) Equation ([5|) gives 

oo 

H(V) = £ U~ 2m+1 J2 c 2m , x V x (5) 

m=l X 

where Vx = ELexK- The coefficients c 2rrti x may be computed from the h 2m ,x in a 
straightforward manner. In particular, they are nonzero only if X is contained in the set 
of sites visited by a 2m step nearest neighbor walk that ends where it begins. The C2 m ,x 
satisfy a bound like (f|). 

For the proof it is useful to split the effective Hamiltonian into three parts. 

H = H, + H 10 + (6) 

H4 contains the terms in eq. ([5]) for m — 1 and 2. H w contains the terms with m = 3,4 
and 5, and contains the terms with m > 5. The fourth order Hamiltonian H4 has 
been computed before 0. In the occupation variables it may be written in the following 
form. 

H 4 = U- x [8 £ V x V y -16Y,Vx] 

<xy>:\x-y\=l x 

+tr 3 [64 e v x v y + iQ J2 v x v y -mJ2 v *+ E c ^v x ] 

<xy>:\x-y\=V2 <xy>:\x-y\=2 % XeE 4 

is the collections of sets X which appear at fourth order and contain a pair of sites 
x, y with \x — y\ = 1. There are only a few such sets and we can compute their coefficients 
explicitly, but as we will see their actual values play no role in the proof. If X e E4 and 
Vx 7^ 0, then there is a nearest neighbor pair x, y with V x V y = 1. This gives a contribution 
of 8U~ l to H4 which is much larger than the order U~ 3 contribution from c^xVx- 

It was shown in || that for configurations with density between 1/5 and 1/4 which 
minimize H4 every 3 by 3 block of sites must look like one of the four configurations 
shown in figure 2 . (We will include a proof of this later.) Watson || made the following 
observation that will play an essential role in our proof. 
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Figure 2 : For densities between 1/5 and 1/4, fourth order perturbation theory 
implies that in a ground state every 3 by 3 square must be one of the above cases, 
up to lattice symmetries. 



Lemma 1 // every 3 by 3 block equals one of those shown in figure 2 (up to a lattice 
symmetry), then the configuration corresponds to a tiling of the lattice by squares and 
diamonds in which every vertex looks like one of the four types shown in figure 3 up to a 
lattice symmetry. 




ABC D 

Figure 3 : For densities between 1/5 and 1/4, every ground state corresponds to 
a tiling of the plane by squares and diamonds in which every vertex must be one 
of the four types shown above. The ions at the center of each of these figures are 
called type A, B, C or D ions, respectively. 



Proof [§§: Note that the 3 by 3 blocks we are using can overlap. Thus the condition that 
every one of these blocks is one of those shown in figure 2 puts many constraints on the 
configuration. We start with a site at which there is an ion. Using the fact that every 3 
by 3 block is one of those shown in figure 2 we work outwards from the initial ion and 
determine all possible configurations in the neighborhood of the inital ion. We find that 
the configuration must look like one of the four cases shown in figure 3 . | 
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It is important to note that the proof of the lemma is local. If we have a configuration 
in which some 3 by 3 blocks are not one of those shown in figure 2 , then we can still 
conclude that in the parts of the lattice where the 3 by 3 blocks are one of those shown in 
the figure the configuration must be given by a square-diamond tiling. We will say that 
an ion is type A, B, C or D if the configuration in the neighborhood of the ion agrees 
with that shown in figure 3 up to a lattice symmetry. Note that in regions where the 
configuration does not correspond to a square-diamond tiling there will be ions that are 
not any of these four types. 

The tenth order Hamiltonian Hiq contains too many terms to list here. However, we 
can organize the Hamiltonian so that we only need the actual values of the coefficients of 
a modest number of them. The values of the coefficients of the terms in H w are important 
only in regions where H 4 is minimized. By the lemma the configuration corresponds to a 
square-diamond tiling in such regions. In these tilings two ions cannot be separated by 
a distance 1, y/2, ^/8 or 3. So we define E w to be the collection of sets that appear at 
tenth order or lower and contain a pair of sites x, y with \x — y\ 2 equal to 1, 2, 8 or 9. For 
X ^ E w the coefficients c 6i x, c 8) x, c 10> x are given in table 1. 



term 


6th 


8th 


10th 


A 


B 


C 


D 


(0,0) 






64 


112 


-704 


1 


1 


1 


1 


(0,0), 


(2,0) 




96 


-1360 


-1440 





1/2 


1/2 


1 


(0,0), 


(2,1) 




216 


-768 


-14000 


2 


2 


2 


2 


(0,0), 


(3,1) 






512 


640 


2 


1 


1 





(0,0), 


(4,0) 






32 


1120 





1/2 


1/2 


2 


(0,0), 


(3,2) 








4000 








1/2 


2 


(0,0), 


(4,1) 








1000 





2 


1 





(0,0), 


(5,0) 








40 


2 





1/2 





(0,0), 


(2,0), 


(1,2) 




-2592 


10560 





1 


1 


2 


(0,0), 


(2,0), 


(4,0) 




192 


1760 











1 


(0,0), 


(2,0), 


(1,3) 






-7040 














(0,0), 


(2,0), 


(3,2) 






-5280 








1 


4 


(0,0), 


(2,1), 


(4,0) 






-3960 





1 


1 


2 


(0,0), 


(2,1), 


(4,1) 






1320 





2 


1 





(0,0), 


(2,1), 


(-1,2) 






-15840 


4 


2 


2 





(0,0), 


(1,2), 


(2,0), (3,2) 






15840 








1/2 


2 



Table 1 : The sixth, eighth and tenth order terms in the Hamiltonian when 
restricted to configurations that correspond to tilings by squares and diamonds. 

Now we are ready to state and prove the main inequalities that will show there is 
phase separation for densities between 1/5 and 1/4 other than 2/9. 
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Theorem 2 Let M be the number of 3 by 3 blocks in which the configuration is not one 
of those shown in figure 2 . Let n A , and no be the number of type A, B, C and D 

ions respectively (figure 3 ). Consider an L by L lattice so there are L 2 sites, and let p 
be the ion density so there are pL 2 ions. Then there are polynomials pi,p2, • • • ,Pe i> n U^ 1 
and functions f\, fz, fz, ft of U~ l such that 

16 

—U~ 3 M <H 4 - p lP L 2 - p 2 L 2 = 0{U~ l )M (7) 

#io - p 3 pL 2 - Pi L 2 = 22AU~ 7 n D + 2UW^n c + 122A0U~ 9 n D + 0(U~ 5 )M (8) 
H w - p 5 pL 2 - p 6 L 2 = 22AU~ 7 n A + 2340C/- 9 n c + \22AW~ 9 n A + 0(f/" 5 )M (9) 
#oo - fipL 2 - f 2 L 2 = 0(U- n )(n c + n D + M) (10) 
#oo - hpL 2 - UL 2 = 0(U- u )(n A + n c + M) (11) 
0{U~ k ) denotes a quantity whose absolute value may be bounded by a constant times U~ k . 

Proof: The inequality for if 4 was proved in 0. We include a short proof for completeness. 
Let B be a 3 by 3 block of sites. (So it contains 9 sites.) Let z be the center of B. Define 

4 16 

H B =U ~\ E V X Vy + U- 3 [l6 Yl V X Vy + — £ V X Vy 

<xy>CB:\x-y\=l <xy>CB:\x~y\=V2 <xy>CB:\x~y\=2 

-lev. -f £ K-f Z 14 + 16+ E C A 

° xeB:\x-z\=l ° x&B:\x-z\=V2 XcB,XeEi " L X 

where mj is the number of translates of X that are contained in B. Now consider Y^b Hb 
where the sum is over all 3 by 3 blocks. (So some of them overlap.) A pair < xy > with 
\x — y\ = 1 is contained in 6 different blocks, a pair with \x — y\ = v2 in 4 different blocks, 
and a pair with \x — y\ = 2 in 3 different blocks. Using these facts we find 

Y,H B = H 4 + (16fT 1 - 64[T 3 )pL 2 + 16fT 3 L 2 (12) 

B 

So to complete the proof we must show that Hb vanishes if B is one of the blocks shown 
in figure 2 and is at least yt^~ 3 otherwise. First note that if the block contains a pair 
of nearest neighbor ions then the second order part of #4 is at least |# _1 - Since the 
fourth order part is smaller by a factor of U~ 2 , this shows that Hb > y#~ 3 f° r an such 
configurations if U is large enough. Now suppose that the block does not contain any 
nearest neighbor ions. Note that this implies that Vx = for all X £ E4. So we can 
easily compute Hb for such configurations. We find that for all 3 by 3 blocks that do 
not contain a nearest neighbor pair of ions and which are not in figure 2 , Hb is at least 
y£/ -3 . Hb vanishes on the configurations in figure 2 . This completes the proof of ([?]). 
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To prove @ and (^|) we first consider configurations which correspond to a square- 
diamond tiling. If X e E w then Vx = in these configurations. We will consider two 
X's to be equivalent if they are related by a translation, reflection and or rotation. For 
X ^ i?io there are 16 equivalence classes, listed in table 1. For each equivalence class 
we want to write the number of X in the class with Vx = 1 in terms of n A ,n B ,nc,nD. 
Consider the second equivalence class. It contains those X's of the form X = {x, y} with 
\x — y\ = 2. The number of such X with Vx ^ in figure 3 is 0,2,2, or 4 for A,B,C, 
or D, respectively. However, this overcounts the number of X with Vx / 0. Each X is 
counted 4 times. So the number of X's in the second equivalence class with Vx = 1 is 
On a + \n B + |nc + n D . These coefficients 0, |, |, 1 along with the coefficients for all the 
other equivalence classes are given in table 1. We should emphasize that the overcounting 
factor is not always 4. It varies from equivalence class to equivalence class, and in one 
case within the equivalence class. Using table 1 we find that for configurations which 
correspond to a square-diamond tiling, 

H 10 = U~ 5 (A96n A + 5Un B + 5Un c + 592n D ) 

+t/~ 7 (-400n A - 4168n B - 4168n c - 7712n D ) 
+£r 9 (-90704n A - 48664n B - 46324n c + 5616n D ) 



The quantities riA,nB } nc and tld are not all independent. Each ion corresponds to a 
vertex in the square-diamond tiling, so 

n A + n B + n c + n D = pL 2 (13) 

By considering the areas associated with each of the four types of vertices in figure 3 , we 
see that 

9 9 

5n A + ~n B + -no + 4n D = L 2 (14) 

Using these two equations we can eliminate n A and n B from our expression for H\q. The 
result is (H). If we use the two equations to eliminate n B and then the result is @. 

This proves @ and @ for configurations that correspond to a square- diamond tiling. 
Now consider a configuration that does not. The number of sites in the region where it 
does not correspond to such a tiling is at most 9M. Thus the terms in H w that intersect 
this region can contribute at most 0(U~ 5 )M. Outside of this region we can apply the 
above argument. There will be errors at the boundary of the region in which there is 
a square-diamond tiling since there will be incomplete squares and diamonds, but the 
contribution of these errors is also 0(U~ 5 )M. Equations (0) and (|14D are not true for a 



general configuration, but the difference between the right and left side of these equations 
is bounded by a constant times M. This completes the proof of @ and (H). 

Finally, we must prove (HI) and (P~T|) . We only give the proof of (0). A similar 



argument proves (|TT|), o r it may be obtained from (|T0D by showing \n A — ri£> — 2L 2 + 
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9pL 2 \ is bounded by a constant times M. To compare a term c 2m ,xVx in -ffoo with the 
corresponding term for the configurations in figure 1 , we need to know if X is contained 
in a region where the configuration agrees with one of the configurations in figure 1 . So 
we make the following definitions. For a lattice site x and positive integer r, let B r (x) be 
the set of sites such that the I 1 distance from x to y is at most r. Given a configuration, 
let A m be the set of sites with a type A ion such that B m (x) is a subset of configuration A 
in figure 1 . B m is defined analogously using configuration B. Define h 2m ,,x = h 2m ,x/\X\ 
so that 



oo 



#00 = E u ~ 2m+1 E E h 2m , x v x (is) 

m=6 x X-.x&X 

Since V x = if V x = 0, we can restrict the sum over x to sites with an ion. 

Let I denote the sites with an ion. We start by comparing with the sum over a 
subset of the terms in H^: 

oo oo 

l^oo - E u~ 2m+1 E E h 2m , x v x \ = i E u~ 2m+l E E 

m=6 xeA m UB m m=6 xei"\(A m UB m ) X:xGX 

oo 

< J2 U ' 2m+1 \ I \( A m^B m )\c m (16) 

m=6 

using (H]). We need to estimate the size of / \ (A m U B m ). Let x E I \ (A m U 5 m ). Then 
B m (x) is not a subset of configuration A or B in figure 1 . So either (i) B m (x) is not 
a subset of any square-diamond tiling, (ii) B m (x) is a subset of a square- diamond tiling 
which contains a type C or type D vertex, or (iii) B m (x) is a subset of a square- diamond 
tiling which contains both type A and type B vertices. In case (i) B m (x) must intersect a 
3 by 3 block that is not one of those in figure 2 . Note that it is not possible for a tiling to 
contain both type A and type B vertices, but no type C or D, so case (iii) never happens. 
Thus we can associate with each site in / \ (A m U B m ) either a "bad" 3 by 3 block or a 
type C or type D vertex. The number of sites that are associated with the same 3 by 3 
block or type C or type D vertex is bounded by dm 2 for some constant d. Thus 



\I\{A m UB m )\ <dm 2 (n c + n D + M) (17) 



Thus (|16D is 



< E U~ 2m+1 c m dm 2 {n c + n D + M)= 0(U- ll )(n c + n D + M) (18) 

m=6 

Let ca be the energy per site from when the entire finite lattice contains config- 
uration A in figure 1 , and let eg be the same quantity for configuration B in figure 1 . 
Let 

oo 

5e A = E U- 2m+1 a 2m (19) 

m=6 
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If x G A m then the definition of A m implies that 

E h 2m ,xV x = a 2m (20) 

X:xeX 

Thus 

oo oo 

lE^" 2m+1 E E h 2m ,xVx - 5n A e A \ < E U~ 2m+1 E W (21) 

m=6 iei m m=6 x£A\A m 

Since A\A m C/\ (An U B M ), ([L7D implies |A \ A m \ < cm 2 (n c + n D + M). This shows 
that © is 0([/- n )(n c + n D + M). 
A similar argument shows 

00 Q 
\J2 U ~ 2m+1 E E WVx-^e B n B | =0(U- u )(n c + n D + M) (22) 

m=6 xGBm X:x&X 

Combining our bounds on ( |i~6"D and ( 

l^oo - 5e A n A - ^e B n B \ = 0(U' n )(n c + n D + M) (23) 
which proves (|I0|). | 

The main theorem follows easily from the inequalities we have proved by a "varia- 
tional" argument. 

Proof of theorem |l|: First we prove the statement about the ground state for density 
2/9. Configuration B has M = and only type B vertices, so theorem |2| immediately 
implies it is a ground state. Any other ground state must have M = and no type C 
or D vertices. The condition M = implies the configuration is a square diamond tiling. 
The only tiling which contains only A and B vertices and has density 2/9 is configuration 
B in figure 1 . 

Now we turn to the proof of the statements for densities between 1/4 and 1/5. We will 
construct a trial configuration with relatively low energy and then use the inequalities in 
theorem ^| to prove theorem |l[ We give the proof for the case of densities in (1/5, 2/9). 
The proof for densities in (2/9, 1/4) is similar. 

Divide the L by L square into two rectangles and put configuration A in figure 1 on 
one side and configuration B on the other side. The relative areas of the two rectangles 
are chosen to give the desired density. The number of 3 by 3 blocks which do not agree 
with one of those in figure 2 is bounded by a constant times L. Thus if 4 — pipL 2 — p 2 L 2 
is 0{U~ l )L. ric and n D are both zero, so H w — p^pL 2 — p^L 2 is 0(U~ 5 )L and — 
fipL 2 — f 2 L 2 is 0(U~ u )L. Any ground state must have energy no greater than that of 
our trial configuration, so theorem |2| implies that in a ground state 
16 

—U~ 3 M + 22AU~ 7 n D + 2MW^n c + 12240[/ _9 n jD + 0(U~ 5 )M 
3 

+0(U- u )(n c + n D + M)< 0{lJ- l )L 



21) with 



shows 
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which may be rearranged as 
16 

[— U~ 3 + 0(U~ 5 ) + 0(U- n )]M + [224U~ 7 + 12240£T 9 + 0([T u )]n c 
3 

+ [2340fT 9 + 0(U- n )]n c < 0(C/~ 1 )L 

If U is large enough this implies that M + ric + n D is bounded by a constant times U 8 L. 
Now take Ao to be the union of all the 3 by 3 blocks which do not agree with one of the 
configurations in figure 2 together with the union over all type C and type D ions of the 
region associated with the ion shown in figure 3 . In A \ Ao the configuration must be 
a square diamond tiling with no C or D vertices. Note that type A and type B vertices 
cannot be adjacent in the tiling, i.e., separated by a distance \/5. However, this does not 
quite insure that each component contains only type A or type B vertices. For example 
one can have a component that consists of two large regions which are connected only by a 
single line of sites. Then one can have type A vertices on one side of the narrow connection 
and type B on the other side. To eliminate this, redefine Ao to be the original Ao plus all 
sites within a distance d of the original set. If d is chosen large enough we eliminate the 
above problem and every connected component of A \ A will be a square-diamond tiling 
with only type A or type B vertices, | 
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